Surface-to-volume ratio with oscillating gradients 
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Restrictions to diffusion result in the dispersion of the bulk diffusion coefficient. We derive the 
exact universal high-frequency behavior of the diffusion coefficient in terms of the surface-to- volume 
ratio of the restrictions. This frequency dependence can be applied to quantify structure of complex 
samples with NMR using oscillating field gradients and static-gradient CPMG. We also demonstrate 
the inter-relations between different equivalent diffusion metrics, and describe how to calculate the 
effect of restrictions for arbitrary gradient waveforms. 



I. INTRODUCTION 

The universal short-time behavior [1] of the diffusion 
coefficient 



D(t) ~ D 1 



D = D\ t=0 , (1) 



allows one to determine the surface-to- volume ratio S/V 
of restrictions in porous materials [2, 3] and in biological 
tissues [4]. Here Do is the unrestricted diffusion coeffi- 
cient, and d is the effective spatial dimensionality, with 
the factor l/d arising from the orientational average of 
the restrictions assuming their statistically isotropic dis- 
tribution [1]. 

While Eq. (1) has been instrumental in characterizing 
restrictions in a variety of media, the direct measure- 
ment of D(i) with pulse field gradient (PFG) diffusion- 
weighted NMR at millisecond time scales is often techni- 
cally challenging, especially in the in vivo applications. 

A promising way to get into the short-time limit is to 
apply the oscillating gradient (OG) method [5], where the 
diffusion weighting is effectively accumulated over many 
periods of oscillation. In this way, the time scale for the 
diffusivity (the oscillation period) can be much shorter 
than the total acquisition time thus enabling practical 
measurements. A variant of this technique requires a con- 
stant diffusion gradient, where the temporal modulation 
is achieved by applying periodic radiofrequency pulses of 
the CPMG type [5-7]. 

In view of applying the oscillating techniques [5-7] , an 
immediate question is, what exactly should one substi- 
tute for the diffusion time t in Eq. (1)? As t ~ 1/w, where 
oj is the gradient oscillation angular frequency, the right- 
hand side of Eq. (1) must transform in the frequency 
domain to 



D 1 




oo . 
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Quite remarkably, the prcfactor in this expression has 
never been explicitly derived for the OG case. The exist- 
ing analytical results are concerned with a finite number 
of echoes [8-13]. Furthermore, there exists a discrepancy 
between the numerical values of Cd provided by different 
groups [13-17]. 

In this work, we find the prefactor C'd exactly both for 
the OG and CPMG cases [Eqs. (10) and (12) below] in 
the limit of a large number of oscillations. This limit is 
practically applicable for high oscillation frequencies in 
accord with the requirement of short diffusion time for 
the validity of Eq. (1). We show that the exact prefac- 
tor values for the infinite OG and CPMG trains differ by 
less than 1% from each other [Eq. (12)], thereby justi- 
fying the view of the CPMG method as being basically 
equivalent to the OG, and validate the approximate nu- 
merical values found in Rcf. [14] for the CPMG and in 
Ref. [15] for the OG. To derive our result, we utilize the 
recently established equivalence between the PFG and 
OG diffusivitics using the effective-medium description 
of diffusion in disordered materials [18, 19]. 

Here, as in Rcf. [1], we do not take into account the con- 
founding effects of heterogeneous magnetic susceptibility 
or relaxation. These effects generally make the inter- 
pretation of the diffusion-weighted measurements chal- 
lenging [13, 20-22]. For the shortest times, they are less 
relevant; in particular, the effect of the surface relaxation 
at the pore walls can be factored out [13]. However, the 
confounding effects can accumulate over the total acqui- 
sition time of many oscillations, significantly modifying 
the apparent S/V ratio. 



II. METHODS 

In this Section we outline the relations between the 
second cumulant of the diffusion- weighted signal and the 
diffusion characteristics relevant for the PFG and OG 
measurements, valid for any statistically isotropic disor- 
dered medium. In Secton III we will apply these relations 
to the problem in question. 
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V(t) = d?[W(t)] 




ReV(uj) q(t) = 7 / t g(i')dt' 



FIG. 1: General relations between the three diffusion metrics: T>(uj), T>(t) and D(t), and the signal attenuation. 



A. The second cumulant 

We begin with the Gaussian phase approximation [23, 
24] to the diffusion-weighted signal, S, 

- lnS(T) ~ 1 / dhdt 2 q(h) (»(ti)«(ta)) g(f 2 ) , (2) 
1 Jo 

which amounts to keeping the second-order term of the 
cumulant expansion [25]. The signal depends on the to- 
tal duration T of the gradient train g(t), and is a func- 
tional of the diffusion-weighting q(t) = 7 f dt' g(t'), with 
7 the gyromagnetic ratio. The diffusion is characterized 
by the autocorrelation function (vftijvfo)} of molecular 
velocity, an even function of t\ — t% in stationary me- 
dia. As we assumed isotropic diffusion from the out- 
set, v here is the velocity component along the fixed 
direction of the applied gradient. For uniform media, 
(v(t\)v(t2)) = 2DoS(ti — £2), leading to the standard ex- 
pression — lnS = bD with b = J Q q 2 (t)dt. 

Here we will utilize an equivalent, and often more con- 
venient way to represent Eq. (2), in terms of Fourier 

transformed quantities, such as = J Q dte tu>t q(t): 

-bxS(T)~ i f q-w {v-uV u ) q u ■ (3) 

The velocity autocorrelator in the frequency representa- 
tion is defined as = J '_ dr e lUT (v(t + r)u(to)) 
independent of to due to time translation invariance. The 
representation (3) underscores that, knowing the corre- 
lator (v- u v u ), one can evaluate the diffusion- weighted 
signal (2) for any gradient waveform g(t). Conversely, by 
selecting a particular form of q(t) according to its Fourier 



representation q u , one effectively allocates a larger or a 
smaller weight to particular Fourier harmonics (v— u v u ) 
contributing to the measured signal (3). 

There are two advantages of working in the frequency 
representation (3). From the practical standpoint, a sin- 
gle integral in u> is simpler than a double integral in t. 
This reduction is due to the time translation invariance 
not explicitly utilized in Eq. (2). From the fundamental 
standpoint, (v_ UJ v UJ ) is directly related to the dispersive 
diffusivity T>(w) discussed below. 

B. The dispersive diffusivity 

As outlined in detail in Ref. [19], the dispersive diffu- 
sivity T)(lu) is a retarded response function relating the 
temporal Fourier component 3 UiT = — I?(o;)V r v E' tJ . r of the 
current J(t, r) of diffusing particles to that of a lump of 
particle density ^(t, r). This makes T>(u>) a central object 
in the effective medium description of diffusion in disor- 
dered media, as it defines the disorder-averaged diffusion 
equation 

-iu)V u , T = 2?HV^ w , r + O (V^, r ) 

which incorporates the characteristics of the restrictions 
that can be quantified with a bulk measurement. At 
the same time, V{lo) = j^°dte lut V(t) is the Fourier 
transform of the retarded velocity autocorrelator T>(t) = 
6(t) {v(t)v(0)}, with 6(t) a unit step function, cf. Fig. 1 
and Ref. [19]. Therefore, 

{v^ u v u ) = 2ReV(uj). 
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As a result, the knowledge of T>(ui) allows one to find 
the second cumulant contribution to the signal attenua- 
tion for any pulse sequence g(t) via Eq. (3): 



lnS(T)~ [^-q_ u V(uj)q u 



(4) 



Here, only Re2?(w) contributes, as ImD(w), odd in cj, 
yields zero after being integrated with an even function 
l^l 2 . Im2?(oj) does not contain additional information 
as it can be restored using the Kramers - Kronig relations 
[26]. As we show below, it may be useful to work with 
the analytic function T>{uS) rather than with its real part. 

The dispersive diffusivity can be obtained exactly 
from the narrow-pulse PFG diffusion coefficient D(t) = 
(x 2 ) /2t via 



V{uj) = D 



dte lu,t d 2 [tD(t)}, 



(5) 



where D = D(t)\ t=0 (cf. Eq. (D3) in Appendix D of 
Ref. [19]). The three diffusion metrics: the dispersive dif- 
fusivity T>(u); the retarded velocity autocorrelator T>(i); 
and the time-dependent diffusion coefficient D(t) contain 
the same amount of information about restrictions, and 
thus can be expressed via each other [19], as illustrated 
schematically in Fig. 1. 



C. Oscillating gradients 

A comprehensive diffusion-weighted measurement 
must provide a way to obtain the diffusivity T>(u>), or 
the correlator (ti_„« u ), for all oj. From this point of 
view, the OG method, with g(t) = go cos Lu^t, is the eas- 
iest one to interpret, as in the limit of the large number 
N = ujqT/2tt 3> 1 of oscillations, 

™19o , x , v xf , si 

effectively selects the ujq component (v- Uo v Uo ), so that 

(190? T 



In S{T) | 



g(t)=go cos u t 



2u 2 



ReX»(a;o). (6) 



Here we used 5(oj)\ u —o = T/2ir from the Fourier repre- 
sentation of S(lj). As a result, it is Rc2?(a;) that is mea- 
sured via the OG techniques [19]. In the above equation, 
the attenuation over each oscillation period is accumu- 
lated, such that the signal S = cxp ( — b ■ Re T>(u)q)) with 
b = Nbi, b\ = n^go) 2 /ojq. For the dispersive V{uj), the 
b- value alone does not define the measurement: the same 
value, achieved with different oscillation frequencies u>o, 
will yield different results for S. 

Remarkably, the signal S is also sensitive to the initial 
phase if of the oscillation g v (t) = go cos(uiot—f), yielding 



-^S(T)\ gv(t) 



sin 2 ip ■ D(T) 
(7) 



where D(T) ~ L>oo = 27(o;)| u ^.o = D{t)\ t ^oo practically 
is the tortuosity asymptote, since the latter is typically 
reached over the sufficiently long total measurement time 
T. Physically, the initial phase if leads to the admix- 
ture of the PFG attenuation over the time T due to the 
nonzero value of 5 w | w -»o °< sintf, cf. Ref. [27] for ip = ^. 

Equation (6), as well as its more general counterpart 
(7), link the diffusive response function T>{uj) of any 
medium to the OG attenuation with N 3> 1 oscillations. 

The above relations reduce the original problem to 
finding the diffusivity for the system in which the 

PFG diffusion coefficient is given by Eq. (1). It can 
be done either by solving the problem [1] in the u>- 
rcprcsentation, or from a Fourier transform of the re- 
tarded velocity autocorrelator T>(t), or directly from the 
time-dependent diffusion coefficient D(t) such as the one 
in Eq. (1), measured by ideal narrow-pulse PFG, via 
Eq. (5). 



III. RESULTS 
A. Dispersive diffusivity at high frequencies 

Below we find the high frequency limit of T>(ui) corre- 
sponding to Eq. (1) in three different ways, in order to 
demonstrate the inter-relations between the above dif- 
fusion metrics (Fig. 1): (i) directly from Eq. (1) using 
Eq. (5); (ii) from the recent result for the diffusivity re- 
stricted by membranes [18] derived in the frequency rep- 
resentation from the very beginning; and (iii) from the 
velocity autocorrelator T>(t) near a flat impermeable wall 
[28, 29]. 

(i) The most direct way to obtain T>(uj) is to use the 
exact relation (5). In general, this relation allows one 
to find T>(ui) for all w knowing D(t) for all t. However, 
it can be also utilized to relate each term of the expan- 
sion of D(t) for short or for long t to the corresponding 
term of the expansion of for high or for low u>, 

respectively. Such expansions in the (fractional) pow- 
ers of t or 1/t usually have finite (or even zero) conver- 
gence radius, hence their frequency counterparts should 
be used within the corresponding bounds. With all that 
in mind, we substitute the second term of Eq. (1) into 
Eq. (5). Reducing the Fourier integral to the Gamma- 
function by rotating the integration contour tot = e l7T ^ 2 u, 
J^dte^H- 1 / 2 = e *"/4 w -i/2r (|) = e"/ 4 v/VZ;, we find 
the universal high frequency limit of the dispersive diffu- 
sivity 




U) — > 00 , 



(8) 



which directly leads to our main result (10) below. 

Equation (8) is the exact universal high-frequency limit 
of the dispersive diffusivity in the presence of restrictions, 
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valid in the limit in which Eq. (1) applies. Further correc- 
tions in the inverse powers of uj will contain information 
about the permeability and curvature of the barriers, as 
well as the spatial correlations between them [18]. 

(ii) In recent Ref. [18], the problem of diffusion re- 
stricted by flat permeable membranes was considered in 
the frequency representation. At high frequencies, this 
solution is completely equivalent to that of Mitra ct al. 
[1] , as in the latter work the impermeable pore walls arc 
approximated by locally flat randomly oriented planes as 
long as the diffusion length is much smaller than the cur- 
vature radius of the walls. Eq. (8) then follows from the 
lj — > oo limit of D(lo) found in Ref. [18], keeping only the 
O^ 1 / 2 ) term. 

(hi) Finally, in Refs. [28, 29] the correction SV(t) to the 
one-dimensional velocity autocorrelator 9(t) (v(t)v(0)) in 
an impermeable box of size L was expressed as the mean 

1 f\ 

ST>(t) = — / v(t)G(xt 7 x e ;t)v(0) dxodx e dxtdxt+e , 
L Jo 

where e — > 0, G(x2, x±; t) is the exact diffusion prop- 
agator, and the velocity operator v(t) = [(xt+ e — 
Xt)/e}G(x t+e ,x t ;e) -> -2D 5'(x t+e - x t ). Indeed, when 
integrated with any smooth function f(x), v gives 



lim / dxdxo - G(x,Xo;e)f(x) 



2D dxdxo f(x)5'(x - x ) = D o [f(L)-f(0)} 



since G at short times becomes Gaussian, G —> Gq{x — 
x ; e) -> 6(x - xo), with G = -2D d x G . Hence 

for L ^> \fDo~b, where 2/L = S/V. In the last expres- 
sion we used the mirror image result for the propagator 
G(x, xo;t) = Gq(x — xo',t) + Gq(x + xo]t) near a wall. 
Averaging over the orientations in d dimensions leads to 



V{t)=Do9{t) \S(t)- 



1 S s/Dr 



(9) 



Eq. (8) is indeed the Fourier transform of this ex- 
pression (here the first term should be understood as 
lim ?? _i. + o S(t—ri), as explained in Appendix D of Ref. [19]). 
This derivation completes the "triangle" of inter-relations 
in Fig. 1. 



B. Oscillating gradients 

Taking the real part of Eq. (8) we arrive at our main 
result, 



RcX>M ~£> 1-C d 




C d 



1 



dV2 



(10) 



in d dimensions. Note that, for the phase-shifted OG se- 
quence, Eq. (10) defines only the first term in Eq. (7), 
whereas the second term, oc D(T), depends on the 
particular system geometry over large spatial scales ~ 

The above value of Cd contradicts the calculation re- 
ported in Refs. [16, 17], where the corresponding prefac- 
tor 1.11 - — fti 1.41 (presumably for the three-dimensional 
case) is about six times greater than our C3. 

In Ref. [15], the above result for d = 3 was rep- 
resented in the form 4c/(9v / 7r) ■ (S/V)y/DoA c s, where 
A e fi = 1/(4/) = 7r/(2w) was called the effective diffusion 
time [30], and the correction factor c rss 0.73 was evalu- 
ated numerically. The exact value is c = 3/4 according to 
Eq. (10), which indicates a 3% deviation in the numerical 
approximation found in Ref. [15]. 



C. CPMG in a constant gradient 

A closely related measurement technique is the CPMG 
train in the presence of a constant gradient [5-7]. Let 
the interval between successive echoes be 2r, with the 
rf pulses applied at t = r, 3r, 5r, ... (cf. notation of 
Ref. [14]). This is equivalent to the box-shaped oscil- 
lating gradients g(t) alternating between the values ±go 
with the frequency uio = 27r/4r, with 4r being the OG pe- 
riod. The Fourier decomposition of this effective square 
gradient waveform 



7(*) = ^E(- 



k COS UJkt 



Lu k = (2k + l)u 



k=0 



2k + 1 

The corresponding q(t) has the Fourier decomposition 
4750 ^ (~) k 



[6(u + (jJk) - S(u - uj k )} 



Substituting it into Eq. (4), we find the signal accumu- 
lated over a large measurement interval T = 2nr, n> 1: 



lnS(T) 



8(75o) 2 T 



r2,,,2 



OO 

^ (2lP 

fc=0 v 



1)4 



Using the above expression (10), find 



mS(T) = 



Tr 2 ( lgo ) 2 DoT 




where 



Cd = Cd 



£9/2 

S4 



0.99351277 C d . 



(11) 



(12) 



Here s v = Efc ( 2fc + l )~ V = i 1 - ^KM, where C(v) 
is the Riemann (^-function. In particular, S4 = 7r 4 /96. 
The term in the exact result (11) can be written 

-0.186843 Sy/Do^/V. The 



as -CdV2/TrSy/DoT/V 
numerical prefactor here agrees well with the approxi- 
mate numerical limit, —0.19 (Ref. [14]), of the calculation 
for the finite number of pulses performed in the time do- 
main [10, 11]. 
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IV. DISCUSSION 

Our approach shows that the exact prefactor C c i, 
Eq. (10) [and its CPMG modification (12)], is as uni- 
versal and independent on the system geometry, as is the 
corresponding coefficient 4/3e?v / 7r in the original result 
(1). The simplicity and generality of this derivation un- 
derscores the utility of the dispersive diffusivity T>(u>). 
Keeping the diffusivity complex-valued simplifies calcu- 
lations in many contexts [18, 19]; taking its real part or 
relating T>(ui) to D(t) is best left for the very last step. 

We also note that in general, the concept of the effec- 
tive diffusion time for the OG protocols [9, 30] is well 
defined only as an order-of-magnitudc estimate, t ~ 1/ui. 
Indeed, the relations between T>(u>) and D{t) are nonlo- 
cal integral relations in time or in frequency [18, 19], i.e. 
to determine T>{uS) one needs to know D(t) for all t, and 
vice- versa, to determine D(t) one needs to know T>(u>) for 
all lu. For the particular short-time limit (1), it was pos- 
sible to relate the 1/yfHj term in the expansion of T>(cj) to 
the corresponding y/t term in D(t), which may prompt 
one to define some effective diffusion time < e ff = so 
that the relative changes in Eqs. (1) and (10) are the 
same (this happens for /3 = 97r/32). However, this ex- 
act proportionality relation generally does not hold for 
all u, i.e. one cannot define some constant (3 such that 
ReV(uj) = D(t)\ t =p/ U for all t. Instead, one has to use 
the exact integral relations [19] between these quantities. 

The effect of restrictions can be calculated for arbi- 
trary gradient waveform using Eqs. (4) and (10) or their 



time-domain counterparts. In particular, one can use a 
gradient waveform defined as a numerical table in mag- 
net's software. 



V. CONCLUSIONS 

In this work we used the equivalence between the de- 
scription of restricted diffusion in time and frequency do- 
mains, to find the exact high-frequency behavior for the 
frequency-dependent diffusivity in disordered media with 
restrictions, accessible with the oscillating gradient and 
static-gradient CPMG protocols. Our results will allow 
one to determine the surface-to-volume ratio of restric- 
tions using the measurement techniques naturally suit- 
able for the shortest time scales. We also demonstrated 
how the effective medium approach unifies and relates to 
each other different diffusion metrics, such as the velocity 
autocorrelator and the time- and frequency-dependent 
diffusion coefficients. 
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